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ABSTRACT 


A nominally spherical drop is used as a model for a theoretical analysis 
of thennocapillary convection and for estimates of convective flow rates in 
"levitated” melts at zero-g. Since in practice temperature fields and the 
resulting convective flow can be more complicated than the simple vertical 
temperature gradient and the single vortex ring, respectively, the convective 
flow arising from a general steady-state temperature field is analyzed. Ex- 
pressions for the components of a steady velocity field are obtained by 
adapting the analytical method of Miller and Scriven. The vortex rings are 
illustrated by means of typical streamlines for the simpler, more symmetric, 
temperature fields. The circulation time is introduced as a measure of the 
rate of circulation in a convection cell and typical values are given for 
several materials. 

INTRODUCTION 

When buoyant forces are negligible, such as in a space laboratory, 
convective flows may still occur in a liquid as the result of gradients 
in the surface— or inter facial tension. These convective flows are commonly 
called the Marangoni effect although the Marangoni effect includes both this 
phenomenon of convection and the phenomenon of the deformation of a free 
liquid surface (Ref. 1) ... Among the causes of gradients in the surface tension 
are gradients in the concentration and in the temperature along the surface 
of a liquid. 

Concentration and temperature gradients may not be completely eliminated 
in many processes, and in some instances, their presence may be necessary to 
produce ..the convection which is desired. The growth of crystals from a melt 
is a process in which the elimination of convection in the melt is desireable 
because convection produces non-uniform growth conditions, and, thereby, an 
increase in the number of dislocations in the crystal. The evaporative purifi- 
cation of a levitated melt is a process in which the opposite result is 
desired; rapid convection is important here because it increases the rate of 
purification by replenishing the impurity concentration at the stir face and 
because it tends to maintain an uniform composition throughout the melt. 


♦Financial support was provided by NASA under contract W-B ,475 #1 
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The question of how much surface tension driven flows might enhance the 
rate of purification of a levitated melt at zero-g is the motivation for the 
work reported here. As a model which will begin to supply part of the answer 
to this question, we will consider thermocapillary convection— that is, 
convection resulting from a temperature gradient along the surface (Ref. 2) — in 
nearly spherical drops. The temperature field at the surface of the drop 
will be treated in a general way by writing it as an expansion in spherical 
harmonigs. The model will examine one of the convective modes corresponding 
to one of the terms in the expansion of the surface temperature field. The 
equations of motion — the Navi er- Stokes equations — will be solved within the 
assumptions of a steady-state and of creeping flow — terms which are nonlinear 
in the velocity will be ignored. The mathematical analysis will proceed along . 
the lines of the method which Miller and Scriven (3) applied to the oscillations 
of a fluid droplet although here we will not retain the time-dependence of 
their problem. The physical boundary conditions will account for both the 
convective flow and for the deformation of the surface so that both aspects 
of the Marangoni effect will appear in the problem. The velocity field 
components which are obtained from the Miller- Scriven analysis will be used 
to derive the circulation time x whose inverse characterizes the rate of 
circulation within a convection cell. An expression for t will be worked 
out in detail for convection cells having axial symmetry. Estimates of x will 
be given for a variety of materials when the convection pattern is a single, 
axially symmetric vortex ring. Also, relative circulation times will be 
calculated for several higher order, axially symmetric convective modes. 
Illustrations of these convective modes will be given. 

THE TEMPERATURE FIELD 

The temperature field responsible for the convection is considered to be 
a general, non-uniform, but steady, field which can be written as: 

T(r,9,«) « T + S Z £_T® rV (9,*) (o-e,o), (1) 

o n**l m»Q o mn mn 

(2a) 

(2b) 

are spherical harmonics as defined by Morse and Feshbach (4) . The function 
pJJ (cos9 ) is an associated Legendre polynominal . The simple case of a constant 
vertical temperature gradient through the drop — (cos8) « cos 8— is 
a special case of the solution of Young, Goldstein and Block (5) who derived 
the velocity field in a bubble rising in a vertical temperature gradient. 


Y e * P m (cosQ) cos m <p 
mn n 

Y° m P m (cos8) sin m <p 
mn n 
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SOLUTION OF THE NAVIER- STOKES EQUATIONS 

To obtain the velocity field u(r,0,4>) and the hydrodynamic pressure field 
p(r,0,4>) within a drop of radius R, we solve the linearized Navier-Stokes 
equations 2 

uV u * Vp (3a) 

7*u - 0, (3b) 

where y is the coefficient of viscosity. We impose the requirements 

1) of a finite solution at r * 0; 

2) of the kinematic condition u (R) * 0, where u is the radial 

component of the velocity; &d r 

3) of the physical boundary conditions which will be examined in 
the next section. 

Since Eqns. (3a,b) are linear, a general solution can be written as a 
superposition of all the modes. Thus, it is sufficient to find a solution 
for one of the modes (n,m,a) . 

Eqns. <3a,b) can be integrated according to the method of Miller and 
Scriven. This method integrates (3a,b) in terms of u r and the radial component 
of the vorticity, where the vorticity is defined by 

? »7x u. (4) 

The results of the Miller-Scriven method which satisfy requirements M 1 M and 
"2" are 

u (r,8,<f,) (R 2 -r 2 )Y°(6,4>) (5a) 

r mn mn 

C (r,0,$) = B° r n_ V (6,*), (5b) 

r mn mn 

The remaining integration constants A a and B° will be obtained from the 
physical boundary conditions in the next section. 

The velocity components — u Q and u. — can be obtained from u and ? by 
a relation due to Sani (6); * r r 


u * e r u r + [r /n(n+l)3 [V^u^e^xV^c^J < 


where is the radial unit vector, 
v__ = V - e I- 


is the surface gradient operator and R is the operator 


P 3? r • 


(8) 
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Eqn. (4) and the results for u^, u Q , u^ can be used to derived C Q 

The hydrodynamic part p(r,9,4>) of the pressure can be found 
the divergence of (3a) which yields 72p * o. 

ThUS ' p(r,6,4>) * F 0 r*V 0 (0 • 

an mn 

The coefficient P a can also be obtained from (3a) by 
mn 

r* 7 2 u * ? 2 (ru r ) * (r/u) (3p/3r) , 
from which it can be shown that 
F 0 * 2(2n+3)uA a /n. 


and C . . 

<P 

by taking 


(9) 

( 10 ) 


( 11 ) 


PHYSICAL BOUNDARY CONDITIONS 


The balance of stresses at the surface requires: 

1) that the normal stress on the surface due to the hydrostatic pressure 
and to the motion of the fluid is balanced by the surface tension; and 

2) that the shear stresses due to the variation of the surface tension 
are balanced by the fluid motion . 

The deformation of the surface must be included' in these conditions. The 
deformation is assumed to be small, so that 

R + AR * R [1 + e(9,$) ] , (12) 


where AR is the displacement of the surface and 

e<e,$) - E a Y a (8,<j>) (13) 

mn mn 

is a radial strain. The coefficient E 0 also must be obtained from the 

mn 

boundary conditions . 


The normal stress due to the fluid 


-p - p + F 0 rV - 2u(3u /Br)/ 

*rr o mn mn ■ > r r+R 

is balanced by the surface tension produced stress 


(14) 



(15) 


*A simplification of the surface conditions is introduced at this point by 
neglecting the inter facial dilational elasticity and the inter facial shear 
elasticity which contribute an inter facial viscosity term to the normal stress 
equations — see Scriven (7) and Miller and Scriven (5) for discussions of these 
properties . 
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where y is the local value of the surface tension and R^ and R^ are the 
principal radii of curvature* For small deformations of the surface of 
a spherical drop. Lamb (8) has shown that 


+ i i [2 + (n-1) (n+2) E® Y®J , (16) 

R^ R mn mn , 

the surface tension coefficient is expanded in terms of the mean value y q 
and the temperature coefficient y^ as 


1 

*x 


y + y fT 
o T mn 


+ L° E a ]rV , 

mn mn mn 


(17) 


where L° is the coefficient obtained in the expansion of T(RtAR) of the 
defonnecP surface about the temperature T (R) of the undeformed surface and 

T a R n Y a is the (n,m,a)-term in the expansion of the surface temperature field, 
mn mn 

Setting the o hydrostatic pressure p^ = 2y /R, considering terms to 

first order in^Y^ in the normal stress°condition, we obtain the first 

equation for A and E a , 
mn mn 

6uR n+1 A® /n - [Y (n-1) (n+2) + 2Y L a R n ]E® - 2y T°R n . (18) 
mn o T mn sin T mn 

Instead of solving the shear stress conditions directly, it is more 
convenient to take the surface divergence and the surface curl of the force 
on an element of surface. The surface divergence equation is 



v[ l? <*v 


7 u ] 

II r r*=R - 


Eqn. (19) can be simplified since u^J =* 
then, yields a second equation for A a and E a , 

mn mn 


o. 


(19) 

The divergence condition, 


2 (2n+l) yR n+1 A°/n - (n+1) yL% V » (n+l)Y_T® R**.. 

mn T mn mn T mn 

The radial part of the surface curl equation yields the result that b' 


mn 


or C * 0* 
r 

a a 

Eqns. (18) and (20) can be solved simultaneously for A and E , 

■ mn mn 

A° « n(n+l) (n+2)Y Y_T® /2uR D® 
mn o T mn mn 

__C7 o n 

E * y jr r /d , 
mn T mn mn 

where D 0 » (2n+l) (n+2)y +y m L 0 R n . 
mn o T mn 


(20) 


at 0 1 


(21) 

(22) 

( 23 ) 
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THE VELOCITY FIELD AND THE YORTICITY FIELD 

Defining a characteristic rate of flow 

U* * -n(n+l) ( n+ 2)Y rt T«T^R n /2liO^ # (24) 

mn o T ran ran 

the components of the velocity u can be written as 

u (r.9,4) - u" r^a-rV J 0 '*> (25a) 

r mn mn 

u 0 (r,9,<j>) - 0 ^n( ?n " 1/n > [1 “^i )?2] ‘ (3Y L /39) (25b) 

u.(r,9,<j>) - U® (r n-1 /n) U-(^f)r 2 ) (l/sin6) • (3Y a /39) , (25c) 

9 mn n+1 mn 

where r » r/R. 


Although C = 0 within the drop, 
the definition of the vorticity, Eqn. 


C Q and C generally do not vanish. 
(4), and^Eqns. (25a, b r c) , 


Using 


=* z!„ r n (l/sin9) ml/H) 
a mn mn 

* _z !L rn(3Y l/ 30) » 

9 mn am 


where 


mn 


2<U /R) (2n+3)/n(n+3), 
mn 


(26a) 

(26b) 

(27) 


•Hie rate of flow, or speed, is 

v * tu r ^+ u Q+u^] (28) 

where in general, 

n 

u * 2, 2. Zn (n,m,a) (29a) 

r n*l m^O a r 

U 9 S U 9 (n ' m ' a) (29b) 

n 

% * n-l m^O § u *( n ' m ' ff > • ( 29e) 

Thus, the general expression for v can be very complicated. Only expressions 
for the pure modes will be investigated here. 

For the special case of the lowest order mode ( n,m,a ) * (1,0, e) and where 
DqI * 0, we obtain equations for a levitated drop in a constant vertical 
temperature gradient. In this case, the temperature in the drop is 


T * T + T, Rrcos9 
o 1 


( 30 ) 
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which results 


in the 


velocity components 
t(y t * x H/3v) (l-f 2 )cos0 
(Y^R /3y) (l-2f 2 )sin0 
0 . 


Ola) 

(31b) 


These results can be obtained from the solutions of Young, et. al (5) , as 
mentioned in the Introduction. In addition, we find that 



mn 




so that the drop is spherical in the limit ^/Y q 


0. 


02 ) 


SOME CHARACTERISTICS OF AXIALLY SYMMETRIC MOPES 


For an axially symmetric mode (n,0,e) , the expression for the rate of 
flow, Eqn. (28), is 


(r n-1 /n) [n 2 (1-r 2 ) 2 <P_) 2 + (I-~ ? 2 ) 2 <P V ] 1/2 
n n+l a 


(33) 


e 1 

where v « v/U , P is the nth order Legendre polynomial and P is the 

associated Le§indre polynomial of first degree. Eqn. (3 3) can^e used to 

identify the stagnation points since v * 0 at these points. 


Within the drop, n rings of stagnation points about the axis of the drop 
can be readily identified: take r 2 * (n+l/(n+3) and P n * 0; since P has n 
nodes and 0 < $ s 2 tt, n rings have been identified. On the surface n of the 
drop (r * 1) , the rate of flow v * 0 if P 1 ■* 0. The associated Legendre 
polynomial pj has (n+1) nodes, including one at each pole. These two 
stagnation points at the poles and the (n-1) stagnation rings on the surface 
coincide with the hot and cold spots of the temperature field and define 
the boundaries of the convection cells in the drop. If n > 1, v * 0 
when r » 0. That is, flow does not occur through the center of the drop 
when n > 1. Finally, we must consider the possibility of other internal 
stagnation points. That is, are there other points such that 


n 2 (l-* 2 ) 2 <P ) 2 + 
n 


a+3 ~2 2 

lfc sr r 1 




0 ? 


The answer is no. Since each term is either positive or zero, and since 
neither (1-r 2 ) and (1- (n+3/n+l)r 2 ) nor P n and pj vanish at identical values 
of r and cos 8, respectively, there are no other internal stagnation points. 


from 


For an axially symmetric problem, a 


2 , 9 
r sin 


r sin 9 


30 


9ft 

3r 


stream function ^(r ,9) can be obtained 

(34a) 


(34b) 
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where & - u /c£ (a « r,e). For the (n,0,e)-mode, the streamfunctionr- is 

a a ion 

$ m -r n+i (X-r 2 ) sin9P*(x)/n(n+l) , ( 35 ) 

where x » cos6. For the mode (l,0,e), the streamfunction is that for Hill's 
spherical vortex (9) . Stagnation points, indicated by "N", and typical 
streamlines for Hill's spherical vortex are illustrated in Fig. 1 and for 
the modes n * 2, 3, 4, 5 in Figs. 2-5, respectively. The drops are seen in 
vertical cross-sections. The model rings lie in planes perpendicular to 
the axis in each illustration. The straight lines within the circles 
represent the boundaries of the convection cells, and here ^ m 0. The 
positions of the stagnation points are listed in Table 1 for these five 
modes. 


Table 1. Stagnation "Points** for Axially Symmetric Convection Cells in 
Nominally Spherical Drops 


Mode n 

r 

Interior 

9 

Type 


Surface 

9 

(r - 1) 
Type 

1 

1//2 


ring 

0® 


point 





180® 


point 

2 

0 


point 

0® 


point 


/3/5 

54° 44' 8" 

ring 

90® 


ring 



125° 15* 52" 

ring 

180® 


point 

3 

0 

— 

point 

0® 


point 


/2/3 

39° 13 * 54" 

ring 

63® 

26* 6" 

ring 


If 

90° 

ring 

116® 

33* 54" 

ring 


» 

140® 46' 6" 

ring 

180® 


point 

4 

0 

— 

point 

0® 


point 


/5/7 

30® 33' 20" 

ring 

49® 

6* 24" 

ring 


N 

70® 7* 28" 

ring 

90® 


ring 


!» 

109® 52' 32" 

ring 

130® 

53' 36" 

ring 


M 

149® 26* 40" 

ring 

180® 


point 

5 

0 

— 

point 

0® 


point 


1/2/3 

25® V 2" 

ring 

40® 

5' 17" 

ring 


II 

57® 25' 14" 

ring 

73® 

25' 38" 

ring 


I* 

90® 

ring 

106® 

34* 22" 

ring 


t* 

122® 34' 46" 

ring 

139® 

54 * 43" 

ring 


II 

154® 58* 58" 

ring 

180® 


point 
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AXIS 
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Figure 3. Streamlines and Stagnation Points for the (3,0,e)-Mode 
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AXIS 
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AXIS: 
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THE CTRCHTJX'rmw WQEr 

The circulation time t will be introduced in this section as a means of 
characterizing the rate of convective mass transfer from the interior of the 
convection cell to the vicinity of the surface of the drop. The reciprocal 
circulation time t -1 is defined by 

T* 1 » f u*d t/f dS 

* f C-ndSy^T dS, (36) 

where the circulation integral 
£ u*di 

is taken around the boundary of the cell formed by the intersection of the 
cell with a vertical plane containing the axis of the drop, such as any of 
the cell boundaries shown in Figures 1*5. The second equality follows 
from Stokes's theorem, where the surface integrals 


/C*ndS and / dS 


are taken over the region bounded by the circuit of the line integral; n is 
the unit normal vector to the surface and here is equal to e . According 
to Eqn. (36) , the reciprocal circulation time is the average'* 5 vorticity 
of the cell. The circulation time resembles the period of rotation of a 
rigid body. Indeed, if the fluid circulated about* the vortical canter 
as a rigid body, Eqn. {36) would yield 4 tt times the frequency of rotation 
and t would be the period of the rotation reduced by 4ir. Unfortunately, 
such a simple interpretation of thezmocapillary convection cells is not 
possible, but t can still be used to characterize the rate of circulation 
within the cell. 

For the axially symmetric mode (n,0,e), the integrals in Eqn. (36) 
can be performed yielding 

T - (n+2)CP (xM-3? (x’ )]/2Z« (6* -8M (37) 

ns n s n s+x on s+l s 

for the circulation time for the s-cell, where 1 f s 5 n+1, and x # and x* 
are roots of s s+i 

P (x) * 0. (38) 

n 

Cell-s has its vortical center at the stagnation point whose angular 
position 0 g is given by the root x g of 

P (x) * 0? 
n 


(39) 
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0 lies within the range whose lower bound is 0^ and whose upper bound is 0^ + _ . 

T§e circulation time may be either positive or negative, the sign depending 3 
upon the direction of circulation in the cell. Using Eqns. (1) , (24) and (27) , 
we can show that, apart from the geometric factor, the circulation time t depends 
only on the material parameters y and y t and on the temperature gradient ns 
at the surface. 

The circulation times for the mode n » 1 are given in Tables 2-5 for a 
variety of liquids. Here, we assume "unit conditions": Lg^ » 0, and a 
unit temperature gradient, Tg^ » -l°/cm. The assumption of unit conditions 
is indicated for the circulation time by t°, where the superscript "o'* 
designates unit conditions. Table 2 lists circulation times for some liquids 
at room temperature? Table 3, for five liquid metals? Table 4, for four 
molten oxides; and Table 5 for three molten sodium halide salts. 


Table 2. Circulation Times* for Some Liquids at Hoorn Temperature (298.15 K) 

Y i 1 

Material T Viscosity i » 

mN/K.m** mHa.s** N s 


Acetone 

-0.112 (11) 

DC 200, 20 cs + 

-0.062 (5) 

200 os 

-0.065 (5) 

1000 cs 

-0.061 (5) 

Ethanol 

-0.9832 (11) 

n-Hexadecane 

-0.106 (5) 

Krytox 143 az"* - *" 

-0.1 (13) 

Hater 

-0.1477 (11) 


0.316 

(10) 

0.0323 

19 

(5) 

3.5 

193 

(5) 

34 

793 

(5) 

180 

1.092 

(10) 

0.150 

3.086 

(12) 

0.334 


32.4 (13) 4 

0.8904 (13) 


♦temperature gradient at the poles has a magnitude of l°/cm. 
**1 raN/m m i dyn/cm? 1 mPa.s » 1 cp « 10** 2 dyn-s/cm 2 . 

4 Uow Coming silicone oils of the IX 200 series. 

++A perf luoroalkylpolyether . 


0.06908 
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Table 3. 

Circulation Times*- 

-Metals 

Materials 

y t 

Tempr T 


mN/K*m 

K ° 

Aluminum 

-0.356 (10) 

930 



1070 

Copper 

-0.06 (10) 

1370 



1470 

Mercury 

-0.2049 (11) 

290 



470 

Sodium 

-0.09933 (14) 

470 



670 

Tin 

-0.0706 (13) 

500 



870 

♦The temperature gradient at the poles has 

Table 4. 

Circulation Times*- 

-Oxides 

Material 

y t 

Temp. T 


mN/K*m 

X ° 

A1 ° 

<0.1) + 

2400 

4 J 


2600 

B 0 

0.0354 (18) 

1410 


at 720 X 

1670 

GeO 

0.056 

1750 

2 

at 1390 X 

1930 

SiO, 

0.031 

2280 

2 

at 2000 X 

2680 



2820 


♦Temperature gradient at the poles has a 
^Estimated*. 


4.5 (15) 

2.5 


0.14 

0.08 
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Table 5. Circulation Times*- 

Material Y T (18) 

raN/K-m 

-Sodium Halides 

Temp. T 
K ° 

Viscosity (17) 
mPa.s 

|t| ,s 

NaCl 

-0.0719 

1090 

1.38 

0.220 



1150 

1.08 

0.172 

NaBr 

-0.0809 

1060 

1.28 

0.181 



1170 

i.oo 6 

0.142 

Nal 

-0.129 

1030 

1.15 

0.102 



1100 

0.96 

0.085 


♦Temperature gradient at the poles has a magnitude of l°/cm. * 


The steady- state temperature gradient is obtained by balancing the 
conductivie heat flux through the drop against the radiant heat flux away 
from the drop in the cooler hemisphere. Sufficient heat is applied to the 
hotter hemisphere to maintain the steady-state gradient. Since the radiant 
heat flux “T 4 , liquids such as molten metals which have high thermal conduc- 
tivities will only have steep steady-state gradients at high temperatures. 
Although copper has the highest thermal conductivity of the liquid metals 
considered in Table 3, it is the only one of the five metals which is 
considered at temperatures high enough for a rate of radiant heat loss to 
be attained which could produce a steady-state gradient of l°/cm. Thus, 
the circulation times given for copper in the table are the only ones 
which are physically attainable under steady-state conditions. For the 
other metals, the physcially attainable value of t will be longer than 
t?: the attainable steady convective rates will be X less than the rates 
produced by a gradient of l°/cm. 

The long circulation times estimated for the oxides, BO, GeO and 
Si0 2 show the importance of the viscosity in determining whetherthermocapillary 
convection can occur in the liquid. For these three oxides, thermocapillary 
convection is expected to be negligible. However, A1 0 has a lower viscosity 
than these other three oxides and, thus, a shorter circulation time. Since 
temperature gradients much larger than l°/cm can be attained in molten oxides, 
significant thermocapillary convection might be observed in molten alumina. 

The n m 1 mode has been demonstrated on earth by the work of Young, 
Goldstein and Block (5) who observed the balancing of the thermocapillary and 
buoyancy forces on small bubbles in a vertical temperature gradient. The low-g 
of a space laboratory in addition to facilating the levitation of a drop would 
allow the observation of the higher convective modes, in particular, the n * 2 
mode. The circulation could be observed by the movement of a dye or of a 
radiotracer. 
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Relative circulation times for the axially symmetric modes n*2, 3,4,5 are 
given in Table 6. These times are calculated relative to the single vortex 
(n * X) time t£. The sign of the circulation time for a particular cell indicates 
the direction of circulation with respect to that in the single vortex drop. 

Since the convective rate decreases as n increases, the higher modes appear 
to be less desirable than the n * 1 mode for obtaining rapid purification. 

However, the convective rates for some of these higher modes still can be 
large enough for the convective pattern to be studied and for these modes 
to have some utility for the purification of levitated melts in space. 

SOMMAKY 

Equations describing the convective flow velocities in thermocapillary 
convection cells in drops have been derived. When the thermal conditions 
approximate a pure, symmetric mode— that is, when the temperature field can 
be described by a single Legendre polynomial , the streamlines of the vortex 
and the positions of the stagnation points in each cell pan be calculated. 

The circulation time, which is the reciprocal of the average vorticity of 
the cell and which characterizes the rate of circulation, also can be 
calculated for these high symmetry cells. An estimate of the circulation 
time for a liquid drop is an indicator of the probable significance of 
thermocapillary convection in the drop when a temperature gradient is imposed 
across it. 


Table 6. Relative Circulation Times, t° /t£, for Axially Symmetric Modes 
with n > 1 

■* . * 

Mode n Cells s t° /t? 

ns 

2 1 3.968 

2 -3.968 


3 1 9.981 

2 -7.365 

3 9.981 

4 1 19.994 

2 -13.505 

3 13.505 

4 -19.994 

5 1 34.960 

2 -22.691 

3 20.647 

4 -22.691 

5 34.960 


*t° is the circulation time for the s-cell of the nth mode, 
ns 

t® is the circulation time for mode n =* 1? a temperature gradient with a 
magnitude of l°/cm is assumed at the poles. The sign indicates the 
direction of circulation with respect to the n *» 1 mode. 
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